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ON  THE  MEASUREMENT  OF  GEOSTROPHIC  OCEAN 
CURRENTS  BY  (NADIR)  SATELLITE  ALTIMETRY 


1.0  INTRODUCTION 

The  following  discussion  considers  the  problem  of  measuring  geo- 
strophlc  ocean  current  velocities  with  a  single-beam,  pulse- limited, 
nadir-looking  satellite  radar  altimeter  similar  to  those  on  SKYLAB, 

GEOS,  and  SEASAT.  It  is  assumed  that  sufficient  information  is  avail¬ 
able  to  accurately  determine  the  instantaneous,  local,  dynamic  ocean 
heights  as  measured  with  respect  to  the  marine  geoid.  These  measure¬ 
ments  require  precise  knowledge  of  the  local  geoid,  satellite  orbit,  and 
non-current-related  factors  affecting  measurement  of  ocean  surface 
topography.  It  Is  also  further  assumed  that  the  differences  in  satel¬ 
lite  orbit  slope  errors  at  orbit  crossings  are  negligible  in  their 
effects  on  the  derivation  of  the  local  dynamic  ocean  surface  slopes. 

Principal  altimeter  measurement  corrections  include  uncertainties 
of  the  order  of  20  cm  for  both  atmospheric  propagation  and  ionospheric 
effects  [Saastamoinen,  1971;  Davies,  1965;  Thomason,  et  al,  1979]. 
Dynamic  processes  causing  local  deviations  of  the  ocean  surface  topo¬ 
graphy  from  the  marine  geoid  include:  tides  (30  cm  open-ocean,  15  m 
coastal)  [Neumann  and  Pierson,  1966;  Doodson,  1958;  Stoker,  1957]; 
barometric  pressure  gradients  (up  to  1  m);  wind  set-up  and  storm  surges 
(up  to  10  m  in  coastal  regions)  [Neumann  and  Pierson,  1966];  and  density 
gradients  associated  with  current  systems  (1  m  associated  with  strong 
boundary  currents  such  as  the  Gulf  Stream  and  the  Kuroshio  Current) 
[Stommel,  1966;  Stommel  and  Yoshida,  1972].  The  measurement  problem  is 
illustrated  in  Figure  1.1. 

Manuscript  submitted  February  26,  1980. 
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Fig.  1.1  -  Configuration  of  altimeter  measurement 


2.0  MEASUREMENTS/AS SUMPTIONS 

2.1  Assumptions.  For  the  following  discussions  it  is  assumed  that: 

(1)  The  marine  geoid  is  accurately  known. 

(2)  The  satellite  orbit  is  accurately  determined. 

(3)  Data  are  available  for  accurate  determination  of 
ionospheric  and  atmospheric  propagation  corrections 
to  satellite  altimeter  range  measurements. 

(4)  Satellite  radar  altimeter  range  measurement  precision 
has  been  determined. 

(5)  Dynamic  components  of  sea-level  deviations  due  to 
tides,  barometric  pressure,  wind  set-up,  and  storm 
surges  are  accurately  determined  from  models  and 
other  data  input. 

2.2  Measurements.  For  the  following  discussions,  it  is  presumed  that 
the  assumptions  above  are  satisfied  and  that  the  measurement  to  be 
considered  is  the  current-related  dynamic  ocean  surface  topography 
within  the  precision  limit  of  the  radar  altimeter  measurement.  Using 
time  (t)  as  the  parameter  to  define  position  along  the  sub-satellite 
track,  the  current- related  dynamic  ocean  surface  topography  can  be 
expressed  as: 

D(t)  -  o(t)  -  [h(t)-eI(t)-eA(t)]-  [G(t)+T(t)+B(t)+W(t)]  (2.1 

where: 

0(t)  is  the  satellite  orbit, 

h(t)  is  the  radar  range  measurement, 

Gj(t)  is  the  ionospheric  propagation  correction, 
eA(t)  is  the  atmospheric  propagation  correction, 

G(t)  is  the  marine  geoid, 

T(t)  is  the  tide, 

B(t)  is  the  topographic  variation  due  to  barometric 
pressure,  and 

W(t)  is  the  topographic  variation  due  to  wind  set-up 
and  storm  surges. 


and 


Both  0(t)  and  G(t)  are  defined  with  respect  to  the  center  of  mass  of 
Earth.  The  determinations  of:  0(t),  e^Ct),  G(t),  T(t),  B(t), 

W(t)  are  assumed  to  be  absolutely  correct. 

3.0  DETERMINATION  OF  GEOSTROPHIC  OCEAN  CURRENT 

Following  Sverdrup,  et  al.  [1946]  ocean  currents  may  be  divided 
into  three  groups:  (1)  currents  that  are  related  to  the  distribution  of 
density  in  the  sea  (barotropic  and  baroclinic  geostrophic  currents) ,  (2) 
currents  that  are  caused  directly  by  the  stress  that  the  wind  exerts  on 
the  sea  surface  (wind  driven  currents) ,  and  (3)  tidal  currents  and 
currents  associated  with  internal  waves.  In  accordance  with  the  above- 
stated  assumptions,  the  following  discussion  concerns  only  the  geo¬ 
strophic  (density-gradient)  currents. 

The  geostrophic  current  is  determined  by  a  balance  between  the 
Coriolis  force  (acceleration  due  to  the  earth's  rotation  and  the  pres¬ 
sure  gradients  (slope  of  the  isobaric  surfaces)  in  the  ocean.  The 
numerical  value  of  the  ocean-current  velocity  for  a  given  slope  of  an 
isobaric  surface  is  given  by: 


tanS 

V«  -  *  2^5 


(3.1) 


where : 

Vr(^)  *  geostrophic  current  speed  in  m/sec  at  Latitude  <$>, 

g  *  acceleration  of  gravity, 

tanS^  *  slope  of  the  isobaric  surface, 

SR  *  angle  of  the  dynamic  ocean  surface  to  the  local 
horizontal  plane, 

id  *  angular  velocity  of  rotation  of  earth,  and 
t  »  geographic  latitude  which  is  positive  and  negative 
for  northern  and  southern  hemispheres,  respectively. 

In  the  northern  hemisphere  (<b>0) ,  the  current  direction  is  90°  to 
the  left  of  the  direction  of  positive  slope.  In  the  sourthern  hemis¬ 
phere  ($<0),  the  current  direction  is  to  the  right. 
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2  —5 

Taking:  g  ■  9.80  m/sec  ,  and  w  -  7.27x10  rad/sec  in 
Equation  (3.1) 

4  lanS— 

VR«)  -  6.74  x  104  -jjjS-  (3.2) 

At  45°  latitude,  for  example: 

Vr(45°)  ’  9-53  x  104  tanSR  (3.2a) 

and  if  tanSR  *  10~^  (10  cm  height  difference  over  10  km): 

V  (45°)  -  1  m/sec  -  2  kts.  (3.2b) 

K 

If  the  surface  slope  is  expressed  in  dimensions  of  'Sect  the  expres¬ 
sion  for  geostrophic  current  speed  becomes,  with  oj  *  15  sec/sec, 

VR(<j>)  -  0.3268  tanSR/sin<J>  (3.3) 

at  45°  latitude,  for  example, 

Vr(45°)  =  0.462  tanSR  (3.3a) 

and  if  tanSR*=  1  sec 

Vr(45°)  =0.5  m/sec  *  1  kt.  (3.3b) 


4.0  GEOSTROPHIC  CURRENT  INFERRED  FROM  ALTIMETER  MEASUREMENTS  AT 

THE  CROSSING  POINT 

It  is  well  recognized  that  a  single-beam  nadir  radar  altimeter 
cannot  provide  measurement  Of  the  geostrophic  current  vector  from  a 
single  along-track  surface  slope.  The  current  derived  from  a  single 
along-track  surface  slope  is,  at  best,  the  across-track  current  com¬ 
ponent,  let  alone  the  true  geostrophic  current.  The  scheme  discussed  in 
Section  6.0  gives  reasonable  recourse  for  short  term  descriptions  of  the 
geostrophic  current  system. 

Nevertheless,  the  true  local  geostrophic  current  velocities  can  be 
inferred  from  altimeter  measurements  at  the  satellite  subtrack 
crossing  points  with  the  additional  assumption:  the  geostrophic 
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current  system  has  not  experienced  significant  changes  during  the  time 
interval  which  is  required  for  the  satellite  to  complete  these  crossings 
In  the  case  of  SEASAT-A,  it  takes  the  satellite  less  than  11  hours  to 
complete  a  crossing  pair  [Cutting,  et  al,  1977].  Since  most  of  the  time 
scales  for  geostrophic  phenomena  in  the  ocean  are  longer  than  11  hours, 
the  assumption  is  valid,  most  of  the  time.  From  the  two  local  dynamic 
ocean  surface  slopes  derived  from  the  altimeter  measurements  made  at  the 
crossings  the  local  dynamic  ocean  surface  plane  can  be  constructed  from 
which  its  true  slope  vector  can  be  evaluated  for  this  local  plane 
[Sokolnikoff ,1951] .  Then  the  geostrophic  current  velocities  may  be 
determined  from  Equation  (3.1)  with  the  known  directions. 

Let  OA^  and  OA2  be  these  two  slope  vectors  without  separation  at 
origin  0  in  a  Right  Hand  Rectangular  Cartesian  Coordinates  as  shown  in 
Figure  4.1  where 

Ax  -  (tx  cos  Or  -  ip/ 2) ,  t±  sin  Or  -  p/2),  ^  s^. 

A2  »  (t2  cos  Or  +  p/2) ,  t2  sin(rr  +  p/2),  t£  s2), 

and  t^  and  are  the  along-track  horizontal  distance  and  measured 

slope,  respectively.  The  intersection  angle  p  between  the  ascending  and 

the  descending  orbits  is  bisected  by  the  x-axis  which  is,  of  course,  in 
the  local  east-west  direction.  These  intersection  angles  with  respect 
to  latitude  for  selected  orbit  inclinations  are  shown  in  Figure  4.2. 

The  plane  which  contains  points  A^  A2,  and  the  origin  0  has  the 
general  form  of 

ax  +  b  v  +  c  z  +  d  *  0  (4  1't 

OOOO 


with  the  origin  0  in  the  plane  dQ 
z  *  ax  +  by 


0  or 


(4.2) 


INTERSECTION  ANGLE 


Latitude . 

Fig.  4.2  -  Intersection  angle  with  respect  to  latitude  for 
orbit  inclinations  of  90°,  93°,  98°,  108°,  120°  and  135° 
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Substituting  the  coordinates  of  points  and  Aj  into  Equation  (4.2) 

cos(ij//2)  a  -  sin(i|//2)  b  *  -  s^  (4.3] 


cos(ip/2)  a  +  sin(t|//2)  b  ■  -  s^ 
Therefore,  from  Equations  (4.3)  and  (4.4) 

S1  +  s2 

a  "  "  2sin(^/2) 

.  S1  ~  S2 
2sin(4//2) 


(4.4] 


(4.5] 


(4.6: 


With  a  and  b  known.  Equations  (4.2),  the  plane  which  contains  points  A. 
A^,  and  the  origin  0  is 


S1  +  s2 
2cos(^/2)  x 


2sin0i>/~2)  y  +  2  =  0 


(4.7 


The  direction  cosines  of  the  normal  to  this  plane  are 

s1  ■f  s2  s1  —  s2 

2cos 0|//2)  ’  "  2sin(^/2)  ’  1 


and  the  angle,  Y  »  the  normal  to  the  z-axis  is  defined  as 


cosy 


31  +  s2 


2cos(iJ//2) 


S1  "  S2 


|_2sin(^/2  J 


+  1 


-  1/2 


(4.9 


for  acute  angle.  Since  Equation  (4.7)  represents  a  plane,  the  normal 
line  defined  by  (4.8)  can  be  a  family  of  lines.  Let  the  particular 
normal  line  be  chosen  to  be  such  that  it  intersects  z-axis.  The 
line  which  is  perpendicular  to  this  normal  line  and  yet  passes  through 
the  origin  lies  in  the  plane  defined  by  Equation  (4.7)  and,  from  Equa¬ 
tion  (4.9) 


sinYt 


cosy 


r-i +  s2 

2 

4. 

si  ~  s2 

2  1 

+  lj 

2cos(t^/2) 

2s  in  (W  2) 

-  1/2 


where  y  is  the  angle  between  this  line  and  the  z-axis.  The  slope 


(tanS,,)  of  this  line  to  the  xy  plane  can,  thus,  be  calculated  as 
R 


tanS„  =  tany 
k  n 


2cos 


1 2  r  l 

+  S2  S^  -  s2 

is(^/2)J  +  2sin(^/2) 


2  11/2 


9  2  1/2 

(s2  +  s2  ”  2s!  s2  cos^)  /sinij; 


(4.10) 


for  acute  angle.  The  equation  for  this  normal  line,  from  Expression 
(4.8),  is 


X  -  X 

0 

0 

1 

sl  +  82  1  " 

r®i  ■  s2  i 

2cos  0J;/2)J 

[_2sin(4)/2)J 

V  zo)  -  (0, 

0,  2  )  which 
0 

~  X 

y 

'  S1  +  s2  I  " 

Pi  -  S2  1 

2cos(iJ>/2) 

1  2sinf^/2)J 

2-2 


=  2-2 


(4.11) 


Equation  (4.11)  can  be  written  also  as 

'S1  *  S2  1 
2cos(^/2)J 


S1  -  S2 


2sin(<Jj/2) 


x  + 


y=0 


(4.12) 


~31  + 
2cos 


y  + 


\ 

2s: 


-  s 


0JV2) 

1 


(z  -  z0) 


(4.13) 


2sin(^/2) 


(z  -  2  ) 
o 


(4.14) 


Equation  (4.12)  is,  indeed,  the  direction  of  the  line  on  the  xy-plane 
which  is  perpendicular  to  Equation  (4.11)  and  has  a  slope  shown  by 
Equation  (4.10)  with  the  xy-plane.  Therefore,  from  Equation  (4.12), 


s2  '  S1 


s2  +  S1 


cot('^/2) 


1 
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or  the  angle  9  with  the  x-axis  Is 


S  "*  s 

tanQ  =  ■— — ~ — —  cot(4</2)  (4.15) 

s2  +  Sj_ 

The  whole  mathematical  scheme  so  far  can  be  categorized  as: 

(1)  OA^  represents  the  slope  vector  for  the  ascending  orbit  which 
intersects  decending  orbit  at  an  intersecting  ground  angle  of  ip.  The 
decending  orbit  has  the  slope  vector  0A2.  The  reference  system  is  the 
local  Rectangular  Cartesian  Coordinates  with  x-axis  in  East,  y-axis  in 
North  and  z-axis  in  the  direction  of  the  local  gravity.  There  is  no 
separation  between  OA^  and  Z5a2  at  the  origin  0. 

(2)  Equation  (4.7)  represents  the  plane  which  contains  the  slope 
vectors  OA^  and  0A2 . 

(3)  The  true  local  slope  vector  of  the  plane  represented  by  Equa¬ 
tion  (4.7)  has  the  magnitude  tanS,  represented  by  Equation  (4.10),  and 
has  the  true  directional  angle  0  with  the  x-axis  as  represented  by 
Equation  (4.15).  This  angle  6  has  180°  ambiguity.  Nevertheless  it 
should  not  be  difficult  to  figure  out  the  correct  one. 

It  is  essential,  and,  also,  important  to  evaluate  the  effects  of 
the  noises  in  the  measurements  of  s^  and  s2  on  Equations  (4.10)  and 

(4.15).  Let 


Sj^  =  Sx  +  2AS^ 


S2  +  2AS2 


(4.16) 

(4.17) 


where  and  S2  are  the  desired  measurements  and  AS^  and  AS2  are  the 
noises  associated  with  the  measurements  of  and  S2>  respectively.  The 
factor  2  is  for  the  worst  case  which  indicates  the  same  noise  at  the 
origin  with  opposite  sign.  In  space  application  it  wilL  be  further 
assumed  that  the  noise  levels  on  the  measurements  on  each  track  are 
the  same,  i.e., 

ASl  -  AS2  -  AS  (4.18) 

This  is  not  generally  true,  but  serves  to  bound  the  problem. 


11 


And  in  Equation  (4.18),  AS  is  the  Root-Mean- Square  (RMS)  noise  in  the 
measurement.  Therefore,  Equations  (4.16)  and  (4.17)  are 


I 

tanAS 

e 


[ (2AS)2  +  (S1  +  S2)2  -  2(2AS)(S1  +  S2)cosiJ>]1/2 
sin^ 


(4.27) 


tan0R 


S2  -  S1 
S2  +  Sl 


co  tip  , 


(4.28) 


and 


F(S1,  S2,  AS) 


2(2AS) 

S2  +  S2  +  2(2AS) 


(4.29) 


The  terms  tanSD  and  tanS,,  are  the  desired  terms  without  the  effects  of 
R  R  f 

the  measurement  noise.  The  terms  tanAS£,  tanAS£,  and  F(S^,  S2,  AS)  are 
all  introduced  by  non-zero  AS. 


In  view  of  the  nonlinear  characters  in  terms  of  S^,  S2>  and  AS  in 
Equations  (4.23)  and  (4.24)  for  S  and  0  it  is  desirable  to  set  and 
S2  to  zero  in  order  to  demonstrate  the  different  effects  on  S  and  9 
due  to  the  same  AS  with  various  orbit  inclinations.  Thus,  Equations 
(4.25),  (4.26),  (4.27),  (4.28),  and  (4.29)  are  reduced  to 

tanSR  »  0  ,  (4.30) 


tanAS 

e 

I 

tanAS 

e 


2 >J2  AS  (1  -  cos\p)1/2, 
sinij ) 

2AS 

sin^  ’ 


(4.31) 

(4.32) 


tan0R  »  0 


(4.33) 


and 


F (0,  0,  AS)  -  1 

Therefore  Equations  (4.23)  and  (4.24)  become 

AS  (1  -  cosi|01/2 


(4.34) 


tanS  “  tanAS. 


sin^ 
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2AS  sec(4>/2) 


(4.35) 


and 

tan0  -  tan0R  *  0  (4.36) 

Equations  (4.35)  and  (4.36)  show  that  the  local  slope  generated  by  AS 
alone  is  in  either  the  east  or  the  west  direction.  The  conclusions 
presented  by  Equations  (4.35)  and  (4.36)  are  quite  obvious  and  should  be 
used  to  ascertain  the  validities  of  Equations  (4.23)  and  (4.24). 

Nevertheless,  in  a  typical  Western  North  Atlantic  Ocean,  let  * 

-1 'sec'  and  S2  *  -9/10  sec* which  give  a  Gulf-Stream-like  current.  Let 
the  root-mean-square  noise  be  1  sec",  i.e.,  2.5  cm  root-mean-square 
altimeter  precision  with  measurements  separated  by  10  km  [ref.  Figure 
5.1,  Equations  (4.16)  and  (4.17)].  By  using  Equations  (3.1),  (4.23), 
(4.24),  (4.25),  (4.26),  (4.27),  (4.28)  and  (4.29),  for  the  worst  possi¬ 
ble  case  where  AS  assumes  the  same  sign  as  those  of  and  S2  Table  4.1 
can  be  computed  for  orbit  inclinations  93°,  98°  and  108°  where 

ip  «  crossing  angle  between  ascending  and  descending  orbits  as 
defined  in  Figure  4.1  and  given  in  Figure  4.2, 

VR  =  the  real  geostrophic  current  speed  in  m/sec, 

V  **  The  geostrophic  current  speed  in  m/sec  computed  with  the 
non-zero  root-mean-square  noise, 

<$R  -  the  real  direction  of  geostrophic  current  with  respect  to  the 
local  east  in  a  counter  clockwise  direction,  and 

6  *  the  direction  derived  for  the  geostrophic  current  V  in  the 
same  sense  of  direction  as  defined  foa.  6R. 

Figure  4.3  shows  the  plots  of  VR  and  V  with  respect  to  the  latitude  for 
orbit  inclinations  of  93°,  98°,  and  108°  while  Figure  4.4  shows  the 
plots  of  6r  and  6  for  these  VR  and  V.  The  broken  lines  in  Figures  4.3 
and  4.4  represent  V  and  5,  respectively,  which  are  computed  with  the 
Root-Mean- Square  (RMS)  noise  AS,  and  the  solid  lines  represent  V  and 
<SR,  respectively,  which  are,  of  course,  computed  without  the  RMS  noise. 
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However,  several  observations  can  be  made  regarding  Figures  4.3  and 
4.4.  (1)  There  are  nearly  constant  relationships  resulting  from 

the  assumptions  such  as 

VR  -  C  V  (4.37) 

that  can  be  made  between  each  pair  of  VR  and  V  for  various  orbit  in¬ 
clinations.  In  reality,  the  value  of  AS  is  Gaussian-distributed,  and  a 
range  of  C  in  Equation  (4.37)  should  be  evaluated,  for  each  orbit 
inclination,  and  independently  verified  by  selective  ground  truth  pro¬ 
grams.  Nevertheless  orbit  inclination  108°  does  have  the  smallest 
magnitude  for  C  among  the  inclinations  selected.  (2)  The  effects  due  to 
RMS  noise  on  6R  appear  to  bend  SR  toward  local  north  or  south  in 
direction  as  predicted  by  Equation  (4.36).  (3)  The  orbit  inclination 
108°  shows  to  be  the  best  among  the  three  orbit  Inclinations  in  the 
determination  of  reasonable  VR  and  <SR.  Therefore,  this  orbit  incli¬ 
nation  ia  the  best  for  geostrophic  current  determination  for  the  choices 
selected. 

5.0  UNCERTAINTY  IN  SATELLITE  ALTIMETER  MEASUREMENT  OF  OCEAN  SURFACE 

SLOPE 

A  slope  of  10  ^  (10  cm  height  difference  in  10  km)  is  approximately 
equal  to  an  angular  slope  of  2 'sec?'  If  all  of  the  assumptions  of 
Section  2.1  are  satisfied,  the  uncertainty  of  the  altimeter  along-track 
slope  determination  is  limited  by  the  measurement  precision. 

Conservative  expressions  for  along-track  slope  measurement  un¬ 
certainties  are  shown  in  Figure  5.1  for  the  cases  of  5  cm  RMS  and  10  cm 
RMS  altimeter  noise  AS.  The  RMS  slope  uncertainty  is  expressed  in 
dimensions  of'sec^  In  the  estimation,  if  the  RMS  noise  is  5  cm,  and  the 
measurements  are  made  at  2  points  separated  by  10  km,  the  RMS  slope 
uncertainty  is  10  ^ ,  or  2  "secT* which  is  for  the  worst  case  estimation. 

It  has  been  demonstrated  that  proper  filtering  (weighted  averaging) 
of  the  along-track  radar  altimeter  data  can  greatly  reduce  the  along- 
track  slope  measurement  uncertainties.  Uncertainties  as  low  as  0.1'secv 
have  been  obtained  from  SEASAT  altimeter  data  (B.  Chovitz  and  B.  Douglas, 
National  Oceanic  Survey,  N0AA,  private  communications).  The  averaging 
techniques  define  the  space-scale  of  precise  slope  determinations. 
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Fig.  5.1  -  Along* track  slope  measurement  uncertainties 
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Since  the  current-speed  uncertainties  are  directly  proportional  to  slope 
uncertainties,  without  elaborated  discussions  on  the  proper  weighting 
techniques  and  the  spatial  scale  factors,  illustrative  examples  of 
current  speed  determinations  are  based  upon  the  conservative  expressions 
as  shown  in  Figure  5.1  wherever  they  are  appropriate. 

The  uncertainty  of  determination  of  the  components  of  surface  slope 
along-track  may  be  expressed  as  a  product  of  a  latitude-inclination- 
direction-dependent  Error  Factor  with  the  along-track  slope  uncertainty . 
The  latitude-dependence  of  North  and  East  slope  error  factors  for 
selected  orbit  inclinations  of  73®,  100®,  105®,  and  110®  are  shown  in 
Figure  5.2. 

The  slope  component  uncertainties  as  a  function  of  latitude  for  an 
along-track  slope  uncertainty  of  0 . 85''secs  and  an  orbit  inclination  of 
98®  (approximately  sun- synchronous)  are  shown  in  Figure  5.3.  Under 
the  conservative  assumptions  of  Figure  5.1,  a  slope  uncertainty 
of  0. 85''sec" implies  a  distance  separation  of  50  km  for  10  cm  RMS  alti¬ 
meter  noise,  or  a  distance  separation  of  25  km  for  5  cm  RMS  altimeter 
noise. 

6.0  DISCUSSION 

From  Figures  4.3  and  4.4  it  is  seen  that:  (1)  The  uncertainties  in 
the  determination  of  geostrophic  current  rapidly  decrease  in  their 
magnitude  and,  yet,  maintain  the  nearly  constant  relationship  as 
given  by  Equation  (4.37)  for  increasing  values  of  latitude  <j>.  Never¬ 
theless,  a  range  of  values  of  C,  in  Equation  (4.37),  can  be  determined 
by  ground  truth  for  Gaussian-distributed  RMS  noise  AS.  For  each  orbit 
inclination  the  distribution  of  C  and  its  lower  and  upper  bounds  for  the 
values  of  C  can  be  calculated.  With  the  known  C,  AS,  V,  and  6,  the  true 
local  geostrophic  current  parameters  V  and  <5_  can  be  evaluated  by 
Equation  (4.37)  with  Figures  4.3  and  4.4.  (2)  The  effects  due  to  RMS 

noise  on  the  determination  of  geostrophic  current  directions  show  the 
directional  angle  6,  with  respect  to  the  local  east,  to  be  more  in  local 
east-west  oriented  direction  than  that  shown  by  the  true  directional 
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Fig.  5.3  -  Along-track  slope  component  uncertainties  with 
0.85  sec  slope  uncertainty  for  98°  orbit  inclination 
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angle  SR.  (3)  For  **  -1 “secT"  S2  **  -0.9'secJ  and  AS  ■  1  sec  the  best 
orbit  inclination,  among  the  three  studied  for  the  detection  of  geo- 
strophic  current  velocity,  is  108°.  The  second  choice  appears  to  be  98°. 
(4)  It  is  also  important  to  point  out  that,  with  the  assumptions  made 
for  the  derivation  for  geostrophic  current  calculation,  the  results  are 
meaningful  only  for  latitude  20°  and  up. 

As  has  been  demonstrated  from  SEASAT-A  analyses,  an  altimeter 
precision  of  10  cm  is  adequate  to  detect  and  provide  first-order  speed 
estimates  for  Western  Boundary  currents  (e.g.  Gulf  Stream,  Kuroshio) 
which  have  characteristic  speeds  of  the  order  1  m/sec,  and  represent  the 
higher  portion  of  the  range  of  current  speeds  to  be  measured.  The 
precision  goal  for  measurement  of  current  speeds  within  the  National 
Oceanic  Satellite  System  (NOSS)  is  5  cm/sec.  From  Figures  4.3  and  4.4 
and  Equations  (3.3),  (  3.3a),  and  (3.3b),  it  is  seen  that  altimeter 
surface  measurement  uncertainty  must  be  of  the  order  of  O.l^sec^to 
achieve  the  NOSS  current  speed  measurement  goal  for  20°  latitude  and  up. 

The  altimeter  measurement  requirement  for  0. l^sec^precision  is 
equivalent  to  a  slope  precision  of  5x10  ^  (5  cm  in  100  km) .  This  slope 
precision  is  necessary  to  determine  the  speed  of  such  important  features 
as  the  Southern  Ocean  Current  and  the  Laborador  Drift  Current. 

Since  in  reality,  the  assumptions  made  in  Section  2.1  cannot  be  fully 
satisfied,  the  averaging  and  optimizing  techniques  such  as  those  em¬ 
ployed  by  Huang,  et  al  (1978)  should  be  used  so  that  the  effects  due  to 
the  partial  invalidities  of  some  of  the  assumptions  can  be  alleviated. 

For  high  latitudes  a  network  of  in-situ  current  measurements  used 
by  Wunsch  (1977)  can  certainly  be  used  to  reduce  the  uncertainties  in 
the  altimeter-derived  current  measurements.  These  uncertainties 
are  primarily,  caused  by  the  small  crossing  angles  and  the  lack  of 
information  on  the  north  component  of  the  slopes  at  the  high  latitudes. 

The  main  purpose  of  the  in-situ  current  measurement  is  to  provide  the 
boundary  condition  for  the  current  system. 


A  multi-beam  altimeter  (McGoogan  and  Walsh,  1978)  can  provide 
several  off-nadir  measurements  simultanteously,.  The  deployment  of  this 
type  of  instrument  along-track  as  well  as  across-track  in  space  will 
have  several  advantages: 

(1)  The  true  local  ocean  surface  slope  vectors  can  be 
obtained  for  the  individual  track  without  waiting 
for  crossings  to  occur.  Henceforth,  the  true 
local  geostrophic  ocean  current  velocities,  not 
only  the  components  across-track,  can  be  inferred 
from  a  single  pass. 

(2)  .  The  effects  due  to  orbit  altitude  rate  uncertainties 

on  the  ocean  surface  slope  measurements  are  greatly 
reduced  in  comparison  to  those  provided  by  the 
single-beam  nadir  altimeter. 

(3)  The  resolution  of  both  across-track  and  along-track 
measurements  can  be  changed  if  desired. 

Further  airborne  tests  and  analyses  of  this  multi-beam  instrument  are 
required  to  be  pursued  in  order  to  answer  some  of  the  questions  associ¬ 
ated,  always,  with  new  instruments. 

7.0  CONCLUSION 

Orbit  inclination  of  108°  is  the  best  choice  among  the  three 
selected  orbit  inclincations  (93°,  98°  and  108°)  for  the  detection  of 
geostrophic  current  vectors  above  20°  latitude.  The  second  choice  is 
98°  orbit  inclination.  The  relationship  shown  by  Equation  (4.37)  as 

VR-  C  V 

where  is  the  true  local  geostrophic  current  speed,  V  is  the  geo¬ 
strophic  current  speed  evaluated  with  the  RMS  noise,  and  C  is  a  co¬ 
efficient  can  be  used  to  study  the  possible  accuracy  of  recovering  geo¬ 
strophic  current  velocities  from  satellite  radar  altimeter  measurements 


of  varying  precision  and  in  different  orbits. 

The  uncertainty  in  geostrophic  current  velocity  determination  is  a 
function  of  the  altimeter  along-track  slope  measurement  precision. 
Altimeter  slope  measurement  uncertainty  must  be  less  than  0.1  sec  to 
achieve  a  minimum  uncertainty  in  geostrophic  current  speed  determination 
of  5  cm/sec. 
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